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We study primordial nucleosynthesis abundance yields for assumed ranges of cosmological lepton 
numbers, sterile neutrino mass-squared differences and active-sterile vacuum mixing angles. We fix 
the baryon-to-photon ratio at the value derived from the cosmic microwave background (CMB) data 
and then calculate the deviation of the 2 H, 4 He, and 7 Li abundance yields from those expected in 
the zero lepton number(s), no-new-neutrino-physics case. We conclude that high precision (< 5% 
error) measurements of the primordial 2 H abundance from, e.g., QSO absorption line observations 
coupled with high precision (< 1% error) baryon density measurements from the CMB could have 
the power to either: (1) reveal or rule out the existence of a light sterile neutrino if the sign of 
the cosmological lepton number is known; or (2) place strong constraints on lepton numbers, sterile 
neutrino mixing properties and resonance sweep physics. Similar conclusions would hold if the 
primordial 4 He abundance could be determined to better than 10%. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35,+c; 95.30.-k 



I. INTRODUCTION 

Recent developments in observational cosmology may 
allow the primordial abundances of the light elements, 
including deuterium, to become novel probes of the mass 
and mixing properties of light sterile neutrinos. The high 
precision inference of the baryon-to-photon ratio r\ from 
the observed relative acoustic peak amplitudes in the 
cosmic microwave background (CMB) PI H 0] suggests 
a new way to employ Big Bang Nucleosynthesis (BBN) 
considerations. 

Historically, the comparison of the observationally- 
inferred light element abundances with calculated BBN 
abundance yields has been carried out with the intent of 
obtaining the baryon density. Indeed, on account of the 
near-exponential dependence of the deuterium yield on r\ 
in BBN, the observations of isotope-shifted hydrogen ab- 
sorption lines in Lyman limit systems along lines of sight 
to high rcdshift QSO's provide another, independent high 
precision measure of the baryon content of the universe 
0, 0. This value of r\ currently is in good agreement, 
within errors, with the CMB-derived value 0,[3- 

However, these two independent determinations of the 
baryon density depend on new neutrino physics in dif- 
ferent ways. In particular, the BBN deuterium yield de- 
pends, albeit weakly, on the neutron-to-proton ratio and 
the expansion rate at the BBN epoch and these quan- 
tities, in turn, can depend on the mass / mixing properties 
of sterile neutrinos. 

It has been shown recently that the 4 He abundance 
yield in Big Bang Nucleosynthesis (BBN) can be dramat- 
ically sensitive to medium-enhanced active-sterile neu- 
trino flavor transformation in the presence of a signif- 
icant lepton number H ^j. This sensitivity comes 



about through (post-neutrino-decoupling) neutrino fla- 
vor mixing-induced alterations in the v e and v e energy 
spectra. These alterations cause changes in the weak in- 
teraction rates governing the inter-conversion of neutrons 
and protons, and so ultimately they cause changes in 
the neutron-to-proton ratio in BBN. We show here that 
active-sterile neutrino mixing likewise can induce modest 
changes in the deuterium and 7 Li abundance yields. 

This sets up a potentially new avenue for probing or 
constraining the active-sterile neutrino mixing parame- 
ter space: comparison of the value of r\ derived from 
observationally-inferred deuterium on the one hand and 
the CMB-derived value on the other. Though the He 
BBN yield is far more sensitive to alterations of the 
neutron-to-proton ratio than is the 2 H yield, at present 
the prospects for reliable and precise determination of the 
primordial deuterium abundance might be better than 
for helium. 

The primordial helium abundance is likely between 
23% and 26% by mass [U El El- It may be pos- 
sible to do much better than this by adroit attention 
to issues of radiative transfer and compact blue galaxy 
morphology [lil fl5L HtSj . However, as more QSO lines 
of sight become available, the statistics for deuterium 
abundance determinations in quasar absorption line sys- 
tems will improve. Arguably, we may already know the 
primordial deuterium abundance at least as well as we 
know helium^, 0, Ift| . In any case, it is worth explor- 
ing how much leverage deuterium measurements have 
in constraining the parameter space of sterile neutrino 
mass/mixing values and lepton number (s). 

The LSND anomaly is being re-investigated in the 
mini-BooNE experiment |18l Il9| . A positive signal in 
that experiment would indicate active neutrino coher- 
ent flavor transformation at a mass-squared scale sig- 



2 



nificantly different from the atmospheric and solar neu- 
trino mass-squared differences, <5m 2 tm « 3 x 10~ 3 eV 2 
and drriQ « 8x 10~ 5 eV 2 , respectively. Given the Z°- 
width limit, this would immediately imply the existence 
of a light sterile neutrino. 

If this light sterile neutrino (and its helicity-flipped 
partner, or "sterile antineutrino" ) were to completely 
thermalize in the early universe, there could be both 
an increased 4 He yield, which is possibly unwelcome, 
and trouble with CMB- and large scale structure-derived 
bounds on the sterile neutrino rest mass closure contri- 
bution p, s is m m m s m m m mm. invoca- 
tion of a significant net lepton number can suppress the 
production of sterile neutrinos in the epoch when neu- 
trinos scatter frequently (i.e., prior to weak decoupling), 
thereby easing these constraints |29. l30l l31|. However, 
this lepton number will drive active-sterile resonant pro- 
duction of the sterile neutrino (or sterile antineutrino) af- 
ter weak decoupling [9j . Post- weak-decoupling resonant 
sterile neutrino production would leave the active neutri- 
nos and the sterile neutrino with distorted, non-thermal 
energy spectra which can have a significant impact on 
the neutron-to-proton ratio and the 4 He yield 0] . 

To investigate the effects of these spectral distortions 
and resonant transformation scenarios on the 4 He, 2 H 
and 7 Li BBN yields, wc follow the evolution of the neu- 
trino distribution functions in various resonance sweep 
scenarios and self-consistently couple this with a calcula- 
tion of the light element abundances performed with the 
full BBN nuclear reaction network code. (We employ a 
modified version of the Kawano/ Wagoner code described 
in Ref. |32|.) In Section II we briefly outline the physics 
of active-sterile resonance sweep in the early universe. 
We describe our nucleosynthesis calculations in Section 
III. Results are given in Section IV. Discussion and con- 
clusions are given in Section V. 



II. RESONANT ACTIVE-STERILE NEUTRINO 
FLAVOR TRANSFORMATION IN THE EARLY 
UNIVERSE 

Invocation of a significant lepton number as a dodge 
to full population of both helicity states of a sterile neu- 
trino in the early universe 0, |23, H3, |^| will imply at 
least some resonant, medium-enhanced destruction of 
this lepton number and the concomitant production of 
sterile neutrinos 0, EI ■ In Ref. @ this general picture 
of post- weak-decoupling active-sterile resonance sweep in 
the presence of a net lepton number was laid out in de- 
tail. The single channel active-sterile neutrino conversion 
problem posed in Ref. [j| has recently been solved [To| . 

The weak decoupling epoch is where active neutri- 
nos cease to scatter rapidly enough to exchange en- 
ergy effectively with the background plasma. This takes 
place when the temperature is T 3 MeV. Any neu- 
trino energy spectral distortions which develop after this 
epoch will not be entirely erased by scattering and emis- 
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FIG. 1: Example active neutrino distribution function for a 
forced, adiabatic resonance sweep scenario. 



sion/ absorption processes. 

A pre-existing net lepton number in any of the neu- 
trino flavors can drive medium-enhanced active-sterile 
neutrino flavor transformation in the early universe, both 
in the coherent regime after weak decoupling, and in 
the high temperature regime where de-coherence in the 
neutrino field becomes significant [33l |. The lepton num- 
ber residing in the neutrino sector associated with flavor 
a = e, /x, r is defined to be (in analogy to the baryon-to- 
photon ratio rj) 



where n v 



(1) 



and n 7 are the neutrino, antineutrino, 



and photon number densities at some epoch. The po- 
tential lepton number corresponding to active neutrino 
flavor a is 



(2) 



where also /? = e, /i, r. (Note that neither L Va or C a 
are comoving invariants; we quote values of these as- 
suming no dilution from annihilation, i.e., at epoch 
T« 3 MeV.) 

It is convenient to cast the neutrino transformation 
problem in terms of the scaled-energy e = E v /T, instead 
of the neutrino energy E v , because the former quantity is 
a co- mo ving invariant. A Mikheyev-Smirnov-Wolfenstein 
(MSW) |34l l3f| resonance, or neutrino mass level cross- 
ing, will occur between sterile and active neutrinos for 
electron neutrinos or antineutrinos for scaled-energy 



7T 2 Sm 2 cos 26* 
4 V / 2C(3) G F T 4 



1 



(3) 



and for mu and tau neutrinos or antineutrinos for scaled- 
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FIG. 2: Active neutrino distribution function obtained from 
the active-sterile transformation solution in Ref. [Tfj . 
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where Sm 2 is the mass-squared difference appropriate for 
the active-sterile mixing channel, ( (3) « 1.20206 is the 
Riemann-Zeta function of argument 3, Gp is the Fermi 
constant, and Y e — (n e - — n e +)n~ 1 r/ _1 is the net elec- 
tron number per baryon. Neutrinos transform at reso- 
nance if the terms inside the absolute value symbols in 
Eqs. and (@J arc positive; antineutrinos transform if 
these terms are negative. In practice, for the rather large 
lepton numbers we employ, we can neglect the neutrino- 
electron scattering term (second term in the denominator 
within the absolute value) since r] « 6 x 10~ 10 is so small. 

The general picture of MSW resonance sweep (i.e., 
how the resonance scaled-energy e ros depends on 
time/temperature) in the early universe is evident from 
Eqs. J2J and Q). As the universe expands and the tem- 
perature drops, e rcs will increase from zero. If v a neu- 
trinos propagate through the MSW resonance coherently 
and adiabatically, they will be converted to sterile neu- 
trinos with 100% efficiency and, consequently, the lep- 
ton number will be depleted (£ Q will decrease) 0, |33| . 
Ref. showed that the resonance can sweep smoothly, 
continuously, and (most importantly) adiabatically only 
from e rcs = to e ros = e m ax- Here e max is the value of 
the scaled energy where the product eC a is a maximum, 
4 ax * 2C (3) (e £ — +l)£a(e m ax) 0. Ref. @ 
showed that the MSW resonance cannot sweep smoothly 
and adiabatically beyond this point, e rcs = e max . 

If we force the resonance to continue to sweep adi- 
abatically and continuously past e max , completely con- 
verting v a neutrinos in the portion of the v a distribution 
with scaled energy E u /T < e res , we would completely 
deplete the lepton number (C a — 0) when the resonance 
reaches e rcs = e c . . (where "c.o." stands for cut-off). As 



an example, the initial and final distribution functions 
for a Fermi-Dirac v e energy spectrum with degeneracy 
parameter (chemical potential divided by temperature) 
7;, yc = 0.05 and potential lepton number C = 0.368 for- 
tius forced adiabatic sweep scenario is shown in Figure ^ 
The resulting active neutrino spectrum in this case would 
have a low energy "cut". That is, the distribution func- 
tion would be zero for E„/T < e c . Q .. The sterile neutrino 
produced in this scenario v a — > v s would have a distri- 
bution function identical to the original ^.-spectrum for 
E u /T < e c . , but zero for larger values of scaled energy. 

Recently the active-sterile resonance sweep problem for 
£rcs > £max has been solved ^(J- The resonance does 
sweep continuously past e max to (near) lepton number 
depletion, but it does so non-adiabatically. That is, for 
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s is not 100% efficient. The net result 
is that the resonance must sweep to higher energy to 
significantly deplete the lepton number. The resonance 
will sweep adiabatically as before to e ma x, but there will 
be non-adiabatic, incomplete conversion in the v a — > v s 
channel as e rcs sweeps to higher scaled energy, and then a 
return to complete, adiabatic conversion at large values 
of e r cs- This scenario is depicted in Figure [3 for the case 
of i/ e 's with an initial Fermi-Dirac spectrum and lepton 
numbers L Ue = 0.0343, and L v = L Vt = 0.15. 

A distorted, non-thermal v e (or v e ) spectrum will 
change the neutron-to-proton ratio in BBN and hence, 
the light element abundance yields over the case with 
thermal, Fermi-Dirac energy spectra [j| . This is because 
the v e and/or T> e energy spectra determine the rates of 
the neutron and proton inter-conversion processes, 



v e + n ^ p - 
v e + p ^ n - 
n # p + e~ 



(5) 
(6) 
(7) 



For a given initial potential lepton number, removing 
neutrino population at higher scaled energy in the spec- 
trum results in a larger effect on the neutron-to-proton 
ratio. This is because of the significant neutrino energy 
dependence in the cross sections and dsitribution func- 
tions which factor into the rates of the processes in Eqs. 
0, 0, and 0. However, another factor in this behav- 
ior is that more neutrino population lies in the portion 
of the v e or 77 e spectrum for E„/T > e c-0 . for the values 
of potential lepton number of most interest here. 

For a given initial potential lepton number, the actual 
(full solution) resonance sweep scenario will give bigger 
BBN alteration effects than will the artificial forced con- 
tinuous, adiabatic sweep model [lOj . This is shown in 
Figure^ where we plot the fractional change (in percent) 
in the helium and deuterium abundance yields from their 
standard BBN values for the CMB-determined baryon 
density as a function of the mass-squared difference dm 2 , 
characterizing the v e — > v s oscillation channel. We give 
this for both the adiabatic sweep to e c . . scenario and the 
full resonance sweep solution. In these calculations we 
have taken the initial lepton numbers to be L„ c = 0.0343 
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FIG. 3: Percent change in 2 H and 4 He yield from standard 
BBN calculated as a function of 5m 2 for the two active-sterile 
mixing cases. The two cases are forced adiabatic resonance 
sweep to e c . . (Figure and the full solution (Figure 



and L„ fi = L Vt = 0.15. 

In reality we could expect active-active neutrino con- 
version simultaneous with active-sterile transformation. 
This could leave complicated distortion features in the 
energy spectra of all neutrino species 0. It is the spec- 
tral distortions of v e and V e neutrinos which are most 
important for BBN. Active-active neutrino mixing could 
tend to partially fill in the E u /T < e max portion of the v e 
or ~D e spectrum, though this could be offset by continued 
sweep to even higher E v /T. In any case, to be conser- 
vative in our BBN abundance yield estimates we will in 
what follows assume a smooth, adiabatic resonance sweep 
to e c o . for a given potential lepton number. Therefore, 
our calculated abundance yield changes for given lepton 
numbers and sterile neutrino mass/mixing data will be 
(usually) underestimates. 
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FIG. 4: Percent change in primordial 4 He yield relative to 
the standard BBN zero lepton number, no sterile neutrino 
model is shown as a function of potential lepton number 
£ e = 2L Ve + L u + L Vt . In all cases the baryon-to-photon 
ratio is fixed at the CMB-derived value r\ — 6.11 x 10 -10 . 
The central horizontal line is the standard BBN case. The 
light solid line is the case with lepton numbers but with no 
active-sterile neutrino mixing. In this case we set all indi- 
vidual lepton numbers to be equal. The lower dashed line 
is the case with active-sterile neutrino mixing in the forced, 
adiabatic resonance sweep scenario with lepton number distri- 
bution factor = (L Vfl + L Vt )/2L„ c = 1. The lower heavy solid 
line is the same case but in the full resonance sweep scenario. 
The upper dashed line similarly corresponds to the adiabatic 
resonance sweep scenario but with lepton number distribution 
factor 10, while the upper solid line is the same case in the 
full resonance sweep scenario. All active-sterile mixing cases 
here have 8m 2 — 1 eV 2 . The light horizontal band gives the 
23% — 26% range for Y p allowed by observational bounds. 



three n — > p rates; likewise, the three p — > n rates: 

An -^i/ c +n^p+e- "I" ^n+e+ —>p+& e 4" ^n^p+e~ +P C (8) 



III. PRIMORDIAL NUCLEOSYNTHESIS 
CALCULATIONS WITH NEUTRINO SPECTRAL 
DISTORTIONS 

In general, active-sterile resonance sweep will go on si- 
multaneously with the charged current weak interactions 
that set the neutron-to-proton ratio (n/p), as well as the 
strong and electromagnetic nuclear reactions associated 
with BBN. We will have distorted neutrino v e and v e en- 
ergy spectra with distortions that change in time as the 
active-sterile resonance sweeps and active-active neutrino 
flavor transformation proceeds. These new features ne- 
cessitate handling the weak interactions differently than 
in the standard BBN case. 

The standard BBN code, originally written by R. Wag- 
oner j3(J and later revised by L. Kawano [3?ll3^ |. calcu- 
lates the processes in Eqs. |JS}, (JSJ, and (7J by adding the 



Ap Ap_j_ e — — >v e +n 4~ A (/e _|_j ) _ TO _|_ e + + Ap+e _ +y c — >n ■ (9) 

In the standard calculations, the integrands in A„ and 
X p are manipulated and condensed into a shorter, two- 
part integral to save computation time. This requires the 
neutrino and antineutrino distribution functions (as well 
as the electron and positron) be of Fermi-Dirac form. If 
the neutrino degeneracy parameters are zero, the code 
calculates A„ and A p with a series approximation to fur- 
ther cut down computation time. This can lead to an 
erroneous « 0.5% increase in the neutron-to-proton ratio 

mm. 

In order to implement time-dependent neutrino and 
antineutrino spectral distortions, we rewrote the BBN 
code to calculate each weak rate in Eqs. I©, 
and Q independently, with no series approxima- 
tions, and changed the integration variable to neu- 
trino/ antineutrino energy instead of electron/positron 
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FIG. 5: Same as for Fig.ijbut now with 8m 2 = 10 eV 2 . 
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FIG. 6: Same as for Fig. |1] but now with negative lepton 
numbers. 
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FIG. 7: Same as for Fig. |5] but now with negative lepton 
numbers. 



energy. We then modularized the weak rate calculations 
so that any neutrino or antincutrino distribution func- 
tion could be entered. Likewise, our modifications allow 
us to implement any desired time dependence in these 
distribution functions and they also allow us to calcu- 
late consistently the energy density in (and spectra of) 
any sterile neutrinos which are produced. These modi- 
fications are implemented in four modules. One module 
contains the matrix elements and phase space integrands 
for the weak rates. The phase space integrands call a 
second module which contains the v e and v e distribution 
functions. The third module defines the limits of integra- 
tion. Finally, the fourth module calculates the resonance 
energy for active-sterile transformation and calls the in- 
tegrator, which in turn calls the other modules for the 
integrands and integration limits. In addition, the ex- 
pansion rate of the universe at any temperature is calcu- 
lated self-consistently with all active and sterile neutrino 
distribution functions. 

In all of our BBN calculations, we set the baryon-to- 
photon ratio to 77 = 6.11 x 10~ 10 . This corresponds to 
the central value of the CMB acoustic peak amplitude- 
determined WMAP Three Year Mean result, 77 = (6.11 ± 
.22) x 10~ 10 , which, in turn, corresponds to a baryon rest 
mass closure fraction VL^h 2 = 0.0223 ± 0.0008 Q, where 
h is the Hubble parameter in units of 100 km s _1 Mpc -1 . 

In our calculations, for illustrative purposes, we adopt 
neutron lifetime r„ = 887.8 s. The current world av- 
erage [33 for this lifetime is r„ = 885.7 ± 0.8 s, but a 
recent measurement suggests it could be as small 
as r„ = 878.5 ± 0.7 ± 0.3 s. Our adopted r„ is larger 
than all of these and this has the effect of making our 
distorted-neutrino-spectrum calculations underestimates 
of the actual change in the n/p-ratio and, hence, nucle- 
osynthesis yield deviations. Though these differences are 
small, it must be kept in mind that r n remains uncertain 
to a degree. 

With these choices of 77 and r„, our modified version 
of the BBN code calculates the 4 He mass fraction to be 
Y p = 0.2429 and the deuterium abundance relative to 
hydrogen to be D/H = 2.543 x 10~ 5 . Although the cur- 
rent uncertainty in the WMAP-derived baryon density 
is relatively small (~ ±3.6%), it nevertheless translates 
into a ~ ±5.5% uncertainty in the predicted primordial 
value of D/H. This is because the BBN deuterium yield 
is a very sensitive function of r\. As we discuss below, 
the error in D /H stemming from the current error in the 
CMB-determined r\ precludes using the observationally- 
determined deuterium abundance to constrain the sterile 
neutrino physics discussed in this paper. However, the 
higher precision determinations of 77 in the projected Four 
Year WMAP results lead to an uncertainty of ±0.00047 
in £lbh 2 , corresponding to ±0.117 in 77, while the forth- 
coming Planck mission forecasts ±0.00017 in Vi^ti 1 , or 
±0.045 in 77 ^l|. These more precise determinations of 
77 will translate into commensurately better precision in 
the calculated D/H values. As we discuss below, these 
could allow for new neutrino physics and/or lepton num- 
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FIG. 8: Same as Fig. H but now for the percent change in 
D/H, the deuterium abundance relative to hydrogen. Here 
the outer horizontal band gives the standard BBN range in 
D/H corresponding to the uncertainty in baryon-to- photon 
ratio r) for the Three Year WMAP data, while the middle 
and inner bands give the D/H range for the uncertainty in r\ 
for the Four Year WMAP and the projected one year Planck 
result, respectively. The horizontal dashed line shows where 
neutrino spectral distortion plus lepton number will give a 
4 He yield exceeding 26%. Here 8m 2 — leV 2 . 

ber constraints. 

Corrections to the code, such as time-step corrections 
and Coulomb and radiative corrections, have been dis- 
cussed extensively (e.g., Ref. |42^ : Ref. Most cor- 
rections are a small additive factor to the final helium 
abundance Y p and are functions of the chosen time-step, 
r], and t„. Since this work uses set values for these param- 
eters and presents the results in terms of percent change, 
the additive corrections do not contribute to the final 
results. 

It is beyond the scope of this work to present precision 
element abundance predictions. Our goal here is to illus- 
trate the global trends in element production resulting 
from adding in the active-sterile transformation physics. 
Eventually corrections, such as the coulomb correction, 
should be calculated autonomously in the weak rates in 
order to give the < 1% accuracy desired. 

IV. BBN ABUNDANCE YIELDS WITH 
LEPTON NUMBERS AND STERILE 
NEUTRINOS 

Here we describe the results of our calculations of light 
element primordial nucleosynthesis in the presence of sig- 
nificant lepton numbers and active-sterile neutrino flavor 
mixing. The properties of light sterile neutrinos and the 
lepton numbers of the universe could be related [2^, |33| , 
but here we shall vary them independently to gauge ef- 
fects on BBN abundance yields. We therefore have five 
quantities to vary. 
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FIG. 9: Same as for Fig.|Ubut now with 5m' 2 = 10 eV 2 . 

The first of these parameters is the rest mass of the 
sterile neutrino m s or, equivalently, the mass-squared 
difference 5m 2 « m 2 characteristic of active-sterile neu- 
trino flavor mixing. The second quantity is the effective 
2x2 vacuum mixing angle 9 characterizing the unitary 
transformation between an active neutrino, which we will 
take to be electron flavor |i/ e ), and a sterile state \v s ) and 
mass/energy eigenstates \v\) and |^2), 

\v e ) — cos#|^i) + sin#|^ 2 ), 

\v a ) = -sin%i) +cosf5|^ 2 ), (10) 

with corresponding rest-mass eigenvalues m\ and 7712, re- 
spectively, such that 8m 2 = \m 2 — m 2 \. Here we set 
sin 2 29 = 10~ 3 to conform with the LSND results. Be- 
cause the expansion rate of the universe is so slow at the 
epoch of medium-enhanced coherent MSW sterile neu- 
trino production, flavor evolution is likely adiabatic, at 
least for scaled resonance energy e res < e max IsLUcj. As a 
consequence, our nucleosynthesis results will change little 
with variation in 9 so long as sin 2 2f5 > 1CP 5 . 

As outlined above and in Ref. |j, instead of treating 
the complete 4x4 mass/mixing matrix with its many un- 
known mixing parameters, we shall consider the effective 
2x2 conversion channels v e ^ v s or D e ^ v s and adopt 
two different resonance sweep schemes in an attempt to 
bracket the BBN effects of the active-active plus active- 
sterile mixing-induced spectral distortions. To follow res- 
onance sweep we use: (1) continuous, adiabatic sweep to 
lepton depletion at e c . ., and (2) the full solution of Ref. 
jlOj . The latter resonance sweep scheme gives the most 
dramatic alterations in v e or v e energy distribution for a 
given lepton number, but in reality active-active mixing 
Ve. f M , T or v e ^ , as well as v^ r ^ v s or v^ r =f± u s , 
will likely fill in some of the spectral deficits in v e or v e , 
as will post-decoupling inelastic neutrino and antincu- 
trino scattering. By contrast, the continuous, adiabatic 
sweep to e c . . scenario gives conservative underestimates 
of BBN effects @. 
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FIG. 10: Same as for Fig. |S| but now with negative lepton 
numbers. 
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FIG. 11: Same as for Fig. [5] but now with negative lepton 
numbers. 



The remaining three parameters in our BBN calcula- 
tions are the actual lepton numbers. For a — e,fi,T we 
have 



( ^ 



\12C(3) 



where £(3) ~ 1.20206, T v and T 7 are the neutrino and 
plasma temperature, respectively, and the ratio of neu- 
trino chemical potential to neutrino temperature is the 
neutrino degeneracy parameter r\ Va . While rj Va is a co- 
moving invariant, L Ua is not because the ratio T„/T 7 
varies as temperature drops and the entropy initially in 
the seas of electrons and positrons is transferred to pho- 
tons. (After all e ± -pairs have disappeared, (Xi,/T 7 ) 3 = 
4/11.) The lepton numbers given in our figures are 
for T„/T~ — 1. We assume that all neutrinos and an- 
tineutrinos initially have Fermi-Dirac equilibrium energy 
spectra (e.g., the heavy dashed line in Fig. |2J), so that 
Vu a =-Vv a - 

Active-active 3x3 neutrino mixing with the solar and 
atmospheric mass-squared differences has been shown to 
"even up" the lepton numbers of each neutrino flavor 
to within a factor of -10 |H |H E3- This also 
may be true in the 4x4 mixing case. Since we consider 
v e ^ v s or v e ^ D s , the relevant potential lepton number 
is C e = 2L Ve +L Vlx +L Vt . The lepton number distribution 
factor 



factor = 



2U 



(12) 



is an important determinant of spectral distortion. We 
consider values of this factor between 1 and 10. The big- 
ger the value of this factor, the larger will be the spectral 
distortion for a given value of r\ v<! or 7y Pe ,9| . This is be- 
cause the larger the potential lepton number, the larger 
will be, e.g., e max and e c o ., quantities which set the scale 
for spectral distortion. 



A positive v e degeneracy parameter r\ Ve reflects an ex- 
cess of f e 's over P e 's, which will have the effect of driving 
the reaction in Eq. |(5J to the right, thereby lowering n/p 



and the He yield. However, v e 



flavor conversion 



can effectively distort the v & spectrum so as to lower the 
overall v e number density and thereby shift the reaction 
in Eq. (JSJ back to the left, in effect counteracting the v e 
degeneracy and, possibly, leading to an increase in the 
4 He yield over the standard BBN, r) Ve = result. This 
argument can be reproduced in analogous fashion for a 
v e excess (r) Ue < 0) and v e v s flavor conversion. 

Figures 0] and [S] show the percent change in the 4 He 
primordial nucleosynthesis yield relative to the standard 
BBN model with no lepton numbers and no new neu- 
trino physics as a function of (positive) potential lep- 
ton number. In all cases the baryon-to-photon ratio is 
set to the central CMB-derived value rj = 6.11 x 10 -10 . 
With lepton numbers alone, but without sterile neutrino 
mixing and spectral distortion, the neutron-to-proton ra- 
tio is suppressed and the 4 He yield is decreased relative 
to the standard model. This trend is weakened or even 
completely reversed when v e ^ v s generated spectral 
distortion is included in the calculations. In general, the 
spectral distortions generated in the full resonance sweep 
scenario produce bigger increases in abundance yield over 
the lcpton-number-only case than does the forced adi- 
abatic resonance sweep scenario. As discussed above, 
this stems from the tendency in the full resonance sweep 
mechanism to deplete v e population at higher energy in 
the distribution function. From Figures 0] and [S] it is 
clear that the existence of a sterile neutrino could al- 
ter significantly the relationship between predicted 4 He 
abundance and lepton numbers. It is also clear from 
these figures that improvement in the precision of the 
observationally-inferred value of Y p could make for strin- 
gent new constraints on the sterile neutrino parameters 
and lepton numbers. Even with the (likely overly) gen- 
erous range of 23% — 26% for Y p we see that larger ster- 
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ile neutrino masses together with larger positive lepton 
numbers tend to produce too much He. 

Figures and likewise show the deviation in the 4 He 
yield from the standard model value but now for nega- 
tive lepton numbers. These figures are not simply mirror 
images of Figures 0] and [S] This is because negative lep- 
ton numbers will produce distortions in the v e energy 
spectrum. There is a threshold in the v e capture process 
in Eq. 10, while there is no threshold in the v e capture 
process in Eq. (JjJJ. The result is that we must have a dis- 
tortion extending to higher energy in the v e distribution 
than in the v e distribution to produce the same magni- 
tude change in 4 He yield. This trend is obvious in Fig.0] 
where there is no discernible difference between the case 
with lepton number alone and the case with i> e ^ v s 
mixing in the forced adiabatic sweep scenario with lep- 
ton number distribution factor 1. However, in Fig. [7] we 
see that with large enough 8m 2 , negative potential lep- 
ton number, and lepton number distribution factor it is 
possible that the 4 He yield would fall below 23% by mass 
fraction, at least for the full resonance sweep scenario. 

Similar trends are evident in the deuterium yield as 
shown in Figures 03 El E3 an d El These figures are 
essentially similar in overall structure to those for 4 He. 
However, because the deuterium yield is so sensitive to 
baryon-to-photon ratio ry, in these figures we show bands 
of ranges of calculated D/H corresponding to the quoted 
uncertainty ranges in r] for the Three Year WMAP data 
and for the expected r\ uncertainties in the Four Year 
WMAP and the Planck results. 

The general change in D/H relative to standard BBN 
is similar to that for 4 He. In the lepton number only 
cases with no sterile neutrinos a positive potential lepton 
number with its accompanying suppression in n/p results 
in a decrease in the deuterium abundance yield. Again, 
this trend is reversed for large enough spectral distortion. 
In general, bigger increases in D/H are created by larger 
Sm 2 and larger positive lepton numbers plus v e ^ v s 
conversion in the full resonance sweep solution. 

However, given the current uncertainty in r\ it is evi- 
dent from these figures that no meaningful constraints on 
lepton numbers alone or on combinations of lepton num- 
bers and sterile neutrino properties can be obtained from 
measurements of D/H. The <5D/H produced, for example, 
at 8m 2 = 10 eV 2 for C e > 0.15 could exceed the uncer- 
tainty range in deuterium yield stemming from the Three 
Year WMAP uncertainty in r\, but the v e spectral distor- 
tions accompanying this scenario would produce 4 He in 
excess of 26% by mass fraction, i.e., exceeding the obser- 
vational bound (horizontal dashed line). 

By contrast, the considerably smaller uncertainty ex- 
pected in, e.g., the Planck CMB-determined baryon-to- 
photon ratio could allow for rather stringent constraints 
on (or signatures of) either lepton numbers alone or com- 
binations of lepton numbers and sterile and active neu- 
trino mixing parameters. However, as is evident in Fig. 
El realizing statistically significant constraints or signa- 
tures would require that the observationally-derived pri- 
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FIG. 12: Same as for Fig. 0] but now with Sm 2 = 3eV 2 and 
where the dashed curves are for lepton number distribution 
factors as labeled, up to 10. Here the full resonance sweep 
solution is employed. 



8m 2 = 3 eV 2 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
2L V + L, +L V 

v e v u v t 

FIG. 13: Same as Fig. EH but now for deuterium. 

mordial deuterium abundance be known to better than 
5% accuracy. This is problematic as currently there is 
likely a 15% to 30% error in observationally-determincd 
primordial D/H. Prospects for bettering these errors will 
be discussed in the next section. 

Negative lepton numbers produce alterations in D/H 
yield which are qualitatively similar (with reversed 
trends) to those in the positive lepton number regime. 
However, as for 4 He, the presence of the threshold in 
v e + p — > n + e + serves to lessen the overall quantitative 
impact on | (5D/H | of spectral distortions from given val- 
ues of (negative) C e , Sm 2 and lepton number distribution 
factor. This is shown in Figures ITU1 and ITTI 

The light element abundance yield alterations resulting 
from lepton numbers and spectral distortions can depend 
on 8m 2 , L Ua , L v , L Vr , and resonance sweep physics in 
complicated ways. In Figures El an d El we show 5 4 He 
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FIG. 14: Contours of percent change in 4 He yield relative to 
standard BBN for ranges of Sm 2 = — 11 eV 2 and potential 
lepton number — 0.44, all for lepton number distribution 
factor 3 and the full resonance sweep scenario. 



and 5D/H, respectively, for the full resonance sweep so- 
lution and for Sm 2 = 3 eV 2 for many values of the lepton 
number distribution factor [L v + L Vt ) /2L„ e = 1, 1.5, 2, 
2.5, 3, up to 10. We see that there is a fair increase in 
<5 4 Hc and <5D/H with increasing values of this factor un- 
til it approaches 5. In broad brush, this trend comes 
about because larger values of this factor mean relatively 
lower v e degeneracy and resonance sweep to higher v e en- 
ergy (e.g., larger e max ). Both of these consequences tend 
to increase 4 He and D/H yields. 

Likewise, larger Sm 2 values generally imply an earlier 
onset of resonance sweep and spectral distortion devel- 
opment. This, in turn, means a bigger effect on n/p, 
as the rates for the neutron-to-proton interconversion 
processes in Eqs. (J5J), ©, and are faster at earlier 
epochs where temperature and, hence, average lepton 
energies are higher. These trends are evident in Fig- 
ures El an d E3 where contours of <5 4 He and 5D/H, re- 
spectively, are shown as functions of Sm 2 and (positive) 
C e = 2L [/c + L v + L Ur for the full resonance sweep solu- 
tion and for lepton number distribution factor 3. There 
is little dependence of abundance yield deviation on 8m 2 
in either figure at low values of C e . However, for a given 
Sm 2 , increasing C e tends to delay resonance sweep and 
the development of spectral distortions. This can be off- 
set with larger Sm 2 . As a consequence, for larger values 
of C e we see a significant Sm 2 dependence in both <5 4 He 
and <5D/H. 

The results shown in Figures 1141 and 1151 may help to 
indicate where we might expect active neutrino inelas- 
tic scattering to partially erase or modify the spectral 
distortions we calculate in the coherent neutrino propa- 
gation limit. The earlier the onset of resonance sweep, 
the more significant inelastic neutrino scattering will be. 
This regime will generally be where the Sm 2 dependence 
in abundance deviations is weakest, i.e., for parameters 
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FIG. 15: Same as Fig. [T31 but now for the primordial deu- 
terium abundance yield. 
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FIG. 16: Same as Fig. 0] but now for percent change in 7 Li 
abundance. 

in the upper left hand corners of Figures ll4l and ll5l where 
lepton numbers are small and Sm 2 is large. 

Finally, the deviation in 7 Li abundance ( 7 Li/H) yield, 
<5 7 Li, relative to the standard BBN zero lepton num- 
ber, no new neutrino physics case, is shown as a func- 
tion of (positive) potential lepton number C e in Fig- 
ures d^m 2 = leV 2 ) andH3(<5m 2 = 10 eV 2 ) for sev- 
eral cases. The general trends in <5 7 Li are similar to 
those for 5 4 He. In the lepton number only case with 
no sterile neutrinos, increasing C e and, hence, decreas- 
ing n/p suppresses the 7 Li yield. 7 Li is produced at 
this n principally as 7 Be through 3 He(a, 7) 7 Be. How- 
ever, there is a small contribution to 7 Li from direct 
production via 3 H(a, 7) 7 Li and the tritium, 3 H, abun- 
dance during BBN tracks the n/p ratio. Though <5 7 Li 
can be large for several of the cases shown in these fig- 
ures, it is at this time not a good candidate for lep- 
ton number or sterile neutrino constraint. This is be- 
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FIG. 17: Same as Fig. but now for Sm 2 = 10 eV 2 . 

cause there remains considerable controversy surround- 
ing both the observationally-determined primordial 7 Li 
abundance and the astration/production history of 7 Li 
in stars and in the interstellar medium. At the value 
of the baryon-to-photon ratio r\ adopted here, the calcu- 
lated BBN 7 Li/H yield is a factor ~ 3 higher than the 
7 Li/H abundance inferred from the Spite plateau in hot, 
old halo stars |4£| . 

V. DISCUSSION AND CONCLUSIONS 

We have for the first time self consistently and simul- 
taneously coupled the full BBN nuclear reaction net- 
work with medium-enhanced active-sterile neutrino fla- 
vor transformation. One conclusion from our work is 
straightforward: the existence of light sterile neutrinos 
v s which mix with active neutrinos v e , v lll v r could al- 
ter significantly the relationship between primordial lep- 
ton numbers and the BBN light element ( 2 H, 4 He, 7 Li) 
abundance yields. Our work also shows that precision 
predictions of primordial nucleosynthesis yields for given 
neutrino properties and lepton numbers likely will require 
accurate treatment of the evolution of neutrino spectral 
distortion. Depending on neutrino mixing parameters, 
this may require a full 4x4 quantum kinetic equation 
neutrino flavor transformation scheme. 

However, the calculations we have performed here al- 
low us to point out some intriguing trends that eventually 
may provide a means for constraining lepton numbers 
and/or sterile neutrino properties. In particular, though 
we find that positive (negative) potential lepton num- 
ber causes the BBN 4 He, 2 H, and 7 Li abundance yields 
to be suppressed (increased) relative to the standard 
BBN zero-lepton-number case, these trends are counter- 
acted and even reversed when sterile neutrinos exist and 
medium-enhanced v e ^ v s (p e ^ v s ) takes place. 

An underlying theme of our work is that the BBN 
paradigm, especially as it concerns constraints on neu- 



trino physics, may be changing. We now know to fair 
precision the baryon density from the CMB acoustic 
peak amplitude ratios. The uncertainty in baryon-to- 
photon ratio r\ likely will improve with future observa- 
tions. This trend will culminate in the near term in 
the Planck mission, which is projected to reduce the un- 
certainty in r\ to < 1%. The near elimination of un- 
certainty in this quantity and the neutron lifetime will 
leave the leptonic sector as the principle source of un- 
certainty in conventional, thermal BBN. This is an in- 
triguing development which comes hard on the heels of 
new experimental/observational revelations of neutrino 
mass-squared differences, flavor mixing parameters, and 
CMB- and large scale structure- derived bounds on light 
neutrino mass contributions to closure. 

The active neutrino mixing parameters in particular 
allow 4 He-derived constraints on electron le pto n number 
L Ve H3 to be extended to L Vyj and L Vt HlEl|4|. Cer- 
tainly, improvements in precision in the observationally- 
inferred helium abundance Y p will translate directly into 
improved constraints on lepton number |27| . This is ob- 
vious in Figures El for the value of C e = 4L Uc at 
which the "Lepton Number Only Case" line crosses the 
limits of the horizontal band ( 4 He mass fraction 23%). 
This occurs for L Ve = £ e /4 w 0.045, corresponding to 
v e degeneracy parameter r\ Vis sa 0.07. This is, of course, 
a crude upper limit, corresponding to our adopted range 
for the primordial helium abundance. 

If the uncertainty in the observationally-inferred deu- 
terium abundance D/H could be improved significantly, 
it could become competitive with 4 He as a probe of lep- 
ton number and new neutrino physics. This is not the 
case currently. The statistical uncertainty in D/H is be- 
tween 15% and 30%, as derived from the isotope-shifted 
hydrogen absorption lines observed in Lyman Limit and 
Damped Lyman Alpha absorption systems along lines of 
sight to high redshift QSO's |3, II3- As Figures HiHTDJ 
and llll show. we would need to get the error in D/H down 
to < 5% to enable deuterium to provide constraints on, 
or signatures of, e.g., sterile neutrinos. 

This is problematic in the short term, but not an in- 
conceivable eventuality in the longer term. The current 
deuterium abundance stems from of order a half dozen 
QSO absorption systems. With the projected increase 
in the number of 8m class telescopes and extensions of 
surveys to the southern hemisphere, we might expect to 
increase the number of "clean" [i.e., no interloper cloud) 
Lyman-a absorption systems by a factor of 2 or so. This 
likely will not be good enough. However, the advent of 
30m class telescopes could give many more QSO lines 
of sight and this might provide for much higher preci- 
sion deuterium abundance determinations. This could 
be useful for probes of primordial baryon inhomogeneity 
or the star formation and chemical evolution history of 
the early universe. However, in this paper we point out 
that it also could be useful for constraints on new neu- 
trino physics, especially as regards lepton numbers plus 
sterile neutrinos. 
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If, for example, we knew a priori the sign of the pri- 
mordial lepton number, then there may be a signature for 
sterile neutrinos or at least a means for better constrain- 
ing their properties. A positive lepton number would be 
expected to give a suppression of 2 H and 4 He relative to 
the zero-lepton-number, standard BBN predictions based 
on the CMB-derived baryon density. However, a light 
sterile neutrino could reverse this and give an increase 
in these abundances over the standard case. In broad 
brush, it works in an opposite sense for a negative lepton 
number, though the effect is less dramatic on account of 
the threshold in v e + p — > n + e + . 

Likewise, if experimental evidence exists for active- 
sterile neutrino mixing, then the sign of the deviation 
of 2 H, 4 He, or 7 Li from the standard zero-lepton number 
BBN case may afford a direct measure of the sign and 
magnitude of the potential lepton number. At the least, 
measurements and reasoning along these lines may allow 
for significantly better constraints on lepton numbers in 
this case. 

There is another possibility. Suppose that active- 
sterile mixing properties are measured in the lab and 
precise primordial light-element abundances are obtained 
from observation. Further suppose that the effects 
(abundance deviations from standard BBN) pointed out 
here are not seen. This could indicate that there is 
another mechanism, other than large (l-ki/J > 1CP 3 ) 
lepton number [tj |2!j |^|, operating to suppress ac- 
tive neutrino scattering-induced de-coherence production 



of seas of v s and v s . We know that 
some mechanism must suppress this process because oth- 
erwise there would be a conflict with CMB- and large 
scale structure-derived bounds on light neutrino mass 
contribution to closure (see for example Ref. [9J ). Two 
alternative means for sterile neutrino production suppres- 
sion have been suggested: a low re-heat temperature for 
inflation |54| : and an alteration of neutrino mass/mixing 
properties (i.e., no mass or mixing) at early epochs. 

In any case, the linkage between the light elements and 
new neutrino physics which we have pointed out here 
increases the leverage of new developments in both ob- 
servational cosmology and terrestrial neutrino oscillation 
experiments. 
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